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^^ . Abstract 

We study quasilinear evolutionary partial integro-differential equations of second order 
which include time fractional p-Laplace equations of time order less than one. By means of 
suitable energy estimates and De Giorgi's iteration technique we establish results asserting 
the global boundedness of appropriately defined weak solutions of these problems. We also 
show that a maximum principle is valid for such equations. 
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1 Introduction and main result 



Let T > 0, and i7 be a bounded domain in M.^ . In this paper we are concerned with global a 
^ I priori bounds for weak solutions of quasilinear problems of the form 



d, 



(k*{u — UQ)]—diva{t,x,u,Du) — b{t,x,u,Du), t £ {0,T), x € fl, (1) 



I ^ where Du stands for the gradient of u w.r.t. the spatial variables, k S Lijoc{^+) is a singular 

^— N . kernel, and k * v denotes the convolution on the positive halfline w.r.t. the time variable, that is 

O ; {k*v){t)^ Jl^k{t~T)v{T)dT,t>0. 

We will assume that the kernel k satisfies the following conditions. 

(Kl) k is of type VC, that is (cf. [ll],[in]) k G Lijoc{^+) is nonnegative and nonincreasing, and 
K^ , there exists a kernel I € Lijoc(^+) such that A: * / = 1 in (0, oo). 

c5 ■ (K2) / G Ly{[0,T]) for some g > 1. 

An important example is given by 

fc(i)=gi_„(t)e-^* and /(t) - g,(<)e-^* + /i(l * [5,(-)e-''-])(i), t > 0, (2) 

with a G (0, 1) and /i > 0, see also [23], [lO]- Here gp denotes the Riemann-Liouville kernel 

5/5W=^, i>0, /3>0. (3) 
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In this case, ((T]) amounts to a time fractional equation of order a G (0, 1). Recall that for a 
(sufficiently smooth) function v on IR4., the Riemann-Liouville fractional derivative Dfv of order 
a G (0, 1) is defined by D^v = -^ {gi^a * v). 

Letting p > 1 and fix — (0, T) x 17, we will further assume that the functions a : i^T x R''^"''"'^ — > 
M^ and b : ft-T x R^+^ -^ R are measurable and that they satisfy the structure conditions 

(Ql) (a(i,x,e,r/)|77)> CoM" - co\^\^ ~ Mt,x), 

(Q2) \ait,x,^,ri)\ < Ci^jP-i + cil^r^ + </^i(i,x), 

(Q3) \b{t,x,tT^)\ < C2\vf-^+C2\^\''-' + Mt,^), 

for a. a. {t,x) G Ht, and all ^ G M, 77 G M^. Here Ci,Ci, i — 0, 1,2, are positive constants, and 



(Q4) The parameter 7 lies in the range 

p\ qj ' 



q ^ N 

Nq 



(Q5) The functions ipi, z = 0, 1, 2, defined on fly are nonnegative, ipl ^ G Li{^t)^ and ipo, (^2 S 
Ls{0.t), where 

Here, as usual, q and q' denote conjugate exponents, i.e. - + -^ = 1. 

The function mq — uq{x) is a given data and plays the role of the initial data for the function 
u. We will assume that mq G L2{^). 

Throughout the paper we will further assume that dO. satisfies the property of positive density, 
see Section [2] 

Before describing the main results we give some comments on applications. Problems of 
the form ([T]) arise for example in mathematical physics when describing dynamic processes in 
materials with memory, e.g. in the theory of heat conduction with memory, see |18) and the 
references therein. Time fractional diffusion equations which are obtained by taking k = gi~a in 
([1]) are also used to model anomalous diffusion, see e.g. [16]. In this context, these equations are 
termed subdijfusion equations (the time order a lies in (0, 1)); in the case a G (1, 2), which is not 
considered here, one speaks of superdiffusion equations. We point out that our general setting 
also includes models which describe nonlinear diffusion phenomena. An important special case of 
(dl) is the class of time fractional p-Laplace equations like e.g. ([SJ below. Let us further mention 
that time fractional diffusion equations of time order a G (0, 1) are closely related to a class of 
MontroU- Weiss continuous time random walk models where the waiting time density behaves as 
i-"-i for t -^ cx), see e.g. [H], [13], [S]. 

We say that a function w is a weak solution (subsolution, supersolution) of ([T]) in J7t, 
if u belongs to the space 

Vg,p := {f G L2qi[0,T];L2in))nLpi[0,T];H^{n)) such that 
k*ve C{[0,T];L2{n)), and (fc x= v)\t=o = 0}, 



a(t, X, u, Du) and h{t, x, u, Du) are measurable, and for any nonnegativc test functfon 
with rj\t=T — there holds 



( 



rjt[k * (u — Uo)] + (a(i, x, w, Du)\Drj) — 6(i, x, w, Du)rj\ dx dt — (<, >) 0. (5) 



This definition makes sense, since under conditions (Q1)-(Q5) the integral in ([5|) is finite, by 
Holder's inequality and the parabolic embedding Vq,p '-?> Lr{^T) (see Proposition 12.11 below). 
We point out that ([ij is considered without any boundary conditions, in this sense weak solutions 
of (H)) as defined above are local ones w.r.t. space. We further remark that weak solutions of ([T]) 
in the class Vq^p have been constructed in |24j in the linear case with p = 2. In view of the basic 
energy estimate (see below) and the known results in the case p = 2 the space Vq^p is the natural 
choice for weak solutions in the general case p £ (1, oo). We strongly believe that under stronger 
assumptions on the nonlinearities a and b it is possible to prove the existence of weak solutions of 
([1]) in the class Vq^p by means of the theory of monotone operators and the techniques developed 
in [23]. Notice also that the initial condition u\t=o = uq has to be understood in a weak sense. 
One can show ([24]) that in case of sufficiently smooth functions u and fc * (m — uq), the condition 
(fc * u)\t=o = implies u\t=o = Uq. 

To state our main results we set Tx — (0,T) x d^ and y+ :— maxjy, 0}. By \A\ we denote 
the Lebesgue measure of a measurable set A C M.^ . When we say that a function u e Vq,p 
satisfies u < K a.e. on Tt for some number ii' e R we mean that (u — if )+ G ip([0, T]; i/^ (il)), 
likewise for lower bounds on Tt- This convention allows to formulate our results without extra 
smoothness assumptions on the boundary dVt. Our main result reads as follows. 

Theorem 1.1 Let p > 1, T > Q, and Q, C M^ he a hounded domain. Let the assumptions 
(K1),(K2),(Q1)-(Q5) be satisfied. 

(i) (Suhsolutions) Suppose uq G i2(^) cind that K > is .such that uq < K a.e. in fl. Then 
there exists a constant 

C = C(A^,p,(3',Co,co,C2,C2,7,s,|Z|l,([o^t]),|¥'o + •~P2\L,(nT)^T,\Q.\) (6) 

such that for any weak suhsolution u G Vq,p of ([!]) in VLt satisfying u < K a.e. on Tt there holds 

esssupu < 2(ii: + niax|l,C( / uldxdt)^"^]), (7) 

where r is defined in (Q4) and 

p\^ q) ^ Nq 

(ii) (Supersolutions) Suppose uq G ^2(17) and that K > is such that uq > —K a.e. in D,. Then 
there exists a constant C like in ^ such that for any weak supersolution u G Vq^p of ([I]) in CIt 
satisfying u > —K a.e. on Tt there holds 

essmiu>-2(K + m&xU,C{ I j {-u)ldxdt)^\) . (8) 



Note that ci, Ci, and <y9i, which appear m (Q2), do not play any role in determining the constant 
in (O and dS]), respectively. 

An important special case of ([T]) is the equation 

d^{u - uo) ~ div {\Du\P-^Du) ^ f in Qt, (9) 

with a G (0, 1) and p > 1. We have the following result. 

Theorem 1.2 Let a G (0, 1), p > 1, T > 0, and f2 C K^ be a bounded domain. Suppose that 
uq G -^cxd(^) o,nd that f G Ls{^t) with s > — I — . Let further q > 1 be a fixed number satisfying 

N , N I 
s> — + q > — + -. (10) 

p pa 

Then for any weak solution u G Vq^p of (0) in VLt which is essentially bounded on Vt there holds 

\u\L^(nT) < C{N,p,a,s, |/|L,(f23,),T, |fi|,max{|uo|L^(n),esssup |u|}). (11) 

Tt 

Here the condition on / is sharp, at least in the cases p — 2, a ^ {0, 1) and p > 2, a £ {p' /N, 1), 
as we will show in Section [5j We also remark that the estimate pT|) is stable for a -> 1, that is, 
the constants in the proof remain bounded as a — >■ 1. Hence in this sense we recover well-known 
results for equations like the classical parabolic p-Laplace equation, which can be found in the 
monograph [7]. 

In this paper we further prove that in case of so-called homogenous structures (see Section [6]) 
the weak maximum principle for weak solutions takes the same form as in the classical parabolic 
case. This applies e.g. to equation dH) with / = 0. 

In the literature not much seems to be known concerning a regularity theory for weak solutions 
to ([!]) in the general setting considered in this paper. To our knowledge, the only paper in this 
direction is [20] , where the global boundedness of weak solutions was proved in the case p = 2 
under similar assumptions on the kernel k and the nonlinearities a and b. On the other hand there 
exists a rather well developed regularity theory for degenerate {p > 2) and singular (1 < p < 2) 
parabolic equations of the form ([T]) with dt{k * (u — Uq)) replaced by S^m, see the monograph 
[7] and the references given therein as well as the recent work [5]. This theory includes besides 
local and global Loo -bounds also much deeper results such as Harnack and Holder estimates 
for weak solutions. For the case p = 2 we also refer to [M] and [15]. In the time fractional 
case the situation is much harder due to the nonlocal nature of d" . Recently, a weak Harnack 
inequality was proved for nonnegative weak supersolutions of (|T]) in a special case where p = 2 
and k = gi-a, see [23]. Concerning results in siron^er settings for ([Ij as well as abstract variants 
of it (mostly with p = 2) we refer to [J, [3], [B], [5], [I^, [IB], [21], ^■ 

Our proofs of the global ioo-bounds use De Giorgi's iteration technique and are based on 
suitable truncated energy estimates for weak solutions of ([T]). These estimates are derived by 
combining the techniques from [20] and [7]. A key ingredient is the basic inequality p^ (see 
below) for nonnegative nonincreasing kernels. We further adopt the method of time regularization 
of the equation which goes back to [20] in the weak setting (see also [TO]) and uses the Yosida 
approximations of the operator B defined by Bv — dt{k *v), see Section [2] 

The paper is organized as follows. In Section 2 we collect some preliminary results such as 
the basic inequality ()14p and we explain the time regularization method in more detail. The 



main result is proved in Sections 3 and 4. Section 3 is devoted to the truncated energy estimates 
and in Section 4 we carry out the iteration process. Section 5 gives the proof of Theorem 11.21 
In Section 6 we estabhsh the maximum principle for homogeneous structures, while Section 7 is 
concerned with the case of natural growth conditions. 

2 Preliminaries 

We first discuss an important method of regularizing kernels of type PC. Let k, I S Lijoc{'^+) 
be as in assumption (Kl). For 1 < p < oo, T > 0, and a real Banach space X we consider the 
operator B defined by 



d 
dt 



Bu= -{k*u), D{B) = {ueLp{[0,T];X): k *u £ oH^i[0,T]:X)}, 



where the zero means vanishing a.t t — 0. It is known that this operator is ?Ti-accretive in 
Lp{[0,T];X), cf. [5], [S], [TU]. Its Yosida approximations _B„, defined by i?„ = nB{n + B)~^, n S 
N, enjoy the property that for any u G D{B), one has BnU -^ Bu in Lp([0, T]]X) as n — > cxd. It 
has been shown in [5^ that 

B„w=-(fc„*w), weLp([0,r];X), neN, 
where the kernel fc„ has the representation 

fc„ = fc*/i„, n e N. (12) 

Here /i„ G Xi.;oc(IR+) denotes the resolvent kernel associated with nZ, that is 

K{t)+n{hn*l){t)^nl{t), t > 0, n e N. 

It is further known that the kernels fc„, n G N, are also nonnegative and nonincreasing, and that 
in addition they belong to Hl{[Q,T]), see e.g. [H], [20], [24]. 

We also remark that (Kl) implies that I is completely positive, see e.g. Theorem 2.2 in [3]. 
Consequently, I and hn are nonnegative for all n G N. 

Note further that for any function / G Lp([0,T];X), 1 < r < oo, there holds hn* f ^ f in 
Lp([0, T]; X) as n -^ oo. In fact, setting u = I * f, we have u G D{B), and 



1^^"*")=! 



BnU = — (fc„ * u) = —{k*l*hn*f)^hn*f^Bu^f in Lp{[0, T];X) 



as n ^ oo. In particular, fc„ — >■ fc in Li([0,r]) as n — )> oo. 

We next recall a fundamental identity for integro-differential operators of the form ^[k*u). 
Suppose k G Hl{[Q^T]) and H G C^(M). Then for any sufficiently smooth function u on (0;T) 
one has for a.a. t G (0,T), 

H'{u{t))^^ (k * u)(i) = I (A; * Hiumt) + ( - ff(^(t)) + H'iuitJHt)) kit) 

+ I {H{u{t - s)) - H{u{t)) - H'{u{t))[u{t -s)- u{t)]\ [-fc(s)] ds. (13) 



This follows from a straightforward computation, see also [30], (33]. An integrated version of 
pS]) can be found in [TTJ Lemma 18.4.1]. 

Equation (J13p is highly important for deriving a priori estimates for problems of the form ([1]) . 
In this paper (see also HO]) we will apply it to the functions H+{y) = |(y+)^ and H_{y) = |(y-)^ 
defined for y € M. Here y_ := min{y,0}. Evidently, H± G C^(IR) with derivative H'^{y) = y±, 
y £ K. If the kernel k belongs to Hl{[0,T]) and is nonnegative and nonincreasing, then it follows 
from (|T3)) and the convexity of H± that for any function u e L2{[0,T]), 

u{t)±l{k*u){t)>^^^[k*{u±f){t), a.a.te(0,r). (14) 

The following two lemmas concerning the geometric convergence of sequences of numbers will 
be needed for the De Giorgi iteration arguments below. The first one is contained, e.g., in [T31 
Chapter II, Lemma 5.6], see also jT] Chapter I, Lemma 4.1]. Its proof is by induction. 

Lemma 2.1 Let {Yn}, n = 0, 1, 2, . . ., be a sequence of positive numbers, satisfying the recursion 
inequality 

K+i < C&"r„i+", n = 0,1,2,..., 

where C, b > 1 and a > are given numbers. If 

Yo<C-^'°'b-^/"\ 
then 

in particular y„ — )■ as n -^ oo. 

Lemma 2.2 Let {Yn}, n = 0, 1, 2, . . ., be a sequence of positive numbers, satisfying the recursion 
inequality 

y„+i < Cb" (r„i+" + Y^+') , n = 0, 1, 2, . . . , 
where C, b > 1 and 6 > a > are given numbers. If 

Yo < (2C)-i/"6-i/"', 

then 

Yn < (2C)-i/"6-i/"' &-"/", n e N, 

and thus y„ — ?> as n ^ oo. 

Proof. The assertion follows directly from the proof of the previous lemma and the trivial estimate 

Yr^+i < 2C6"y„i+", 
which holds whenever Yn < 1, due to the assumption S > a. D 

We conclude this preliminary part with an interpolation result which will be frequently used 
in this paper. 



Let T > and flhe a. bounded domain in M.^ . For p,q > 1 we define the spaces 

Vg,p :=V,,pi[0,T] X r!) = L29([0,T];L2(f)))nip([0,T];ffpi(f])), (15) 

and 

both equipped with the norm 

l"|y,,p([O.T]xn) := l"lL2,([0.T];L2(n)) + I^w|Lp([0,T];H1(O))- 

We wiU assume that dfl satisfies the property of positive density, i.e. there exist (5 G (0, 1) and 
po > such that for any xq G F, any ball B{xq,p) with p < po we have that |il n -6(2:0, p)| < 
S\B{xo,p)\, cf. e.g. [21 Section LI]. 

Proposition 2.1 There exists a constant C = C{N,p, q) such that for every u £ V?p([0, T] x Vt) 
there holds 

T 

u(t, x)\^ dx dt 



Jn 

< 

Jn ' ^Jo ' Jn 



1- i + ^ 1 



C"^(/ \Du{t,x)\Pdxdty U { \u{t,x)\^dxydt'^ "' , (17) 

w/iere 

'^= u]_\)'\ '^^d /3 = _^^e[0,l] 

pV-^ q) ^ Nq ^ q^ N 

Proof. We proceed similar as in [71 Chapter I]). We first consider the case where p > 2N/{N + 2). 
By the Gagliardo-Nirenberg inequality (see e.g. [7, Theorem 2.1]), we have for a. a. t € (0,T) 

k(i,-)U.(o) <Ci(7V,p)|i?«(t,.)li^(a)l"(i,-)li;fa), (19) 

where /? and r are given by (fT8|) . In fact, a short computation shows that ([T8| implies that 



/3 



1 _ 1 

2 r 



J. _ 1 I 1' 

N p ' 2 



which corresponds to condition (2.2) in [7J Theorem 2.1]. Taking the rth power in ([T^. integrat- 
ing over (0,T), and using Holder's inequality yields 

i7/r < r'-|n7;i'3'' h/i^^"'^''' ^ _ (i - PVp 

One verifies that f = 2q, and so (fT7|) is valid. 

If p < 2N/(N + 2), then we have in particular p < N, and thus by [71 Corollary 2.1] 

Ht,-)\L^^^n) < C2iN,p)\Duit,-)\L^^n), a.a. te(0,T). 

N-p 



Using this and the fact that 

.(1-/3) ^^ 



Np- {N - p)rl3 
it fohows by means of Holder's inequality that 



Jn Jo Jn 

< / (/ \u\T^dx) "" ( / \u\^dx) ^ dt 
Jo Jn Jn 

< Ct I ( / l^^r dx)^{( |u|2 dx) ^^^ dt. 

Jo Jn Jn 

As in the first case we may now apply Holder's inequality once more thereby proving (J17p . D 
We remark that r > 2 if and only if p > -^^ ■ 

3 Energy estimates 

The following lemma will be the starting point for all of the a priori estimates derived in this 
paper. It provides an equivalent weak formulation of ([Ij where the kernel k is replaced with the 
more regular kernel fc„ (n G N) defined in (fT2|) . In what follows the kernels /i„, n e N, are as in 
Section H 

Lemma 3.1 Let p > 1, T > and fl C M^ be a bounded domain. Let the assumptions 
(K1),(K2),(Q1)-(Q5) be satisfied and assume that uq G L2(^)- Then u G Vq^p is a weak solution 
(subsolution, supersolution) of (QP if and only if for every nonnegative function ip G Hp(fl) one 
has 

(ipdilkn * {u — uo)] + (^n * a{-,x,u, Du)\Dip^ — [ft,„ * b{-,x,u, l?u)]'0) dx 

= (<, >) a.a. t G (0, T), n G N. (20) 

Proof. The proof is analogous to the proof of Lemma 3.1 in 20 . For the reader's convenience 
we repeat it here. 

We may restrict ourselves to the subsolution case as the remaining cases can be treated 
analogously. 

The 'if part can be seen as follows. Given an arbitrary nonnegative 77 G H2 (^t) satisfying 
ri\t=T — 0, we take in pO)) ip{x) — ri{t, x) for any fixed t G (0, T), integrate from t = to i = T, 
and integrate by parts w.r.t. the time variable. Sending then n -^ 00 yields ([S]); here we use the 
approximating properties of the kernels /i„ described in Section [21 

To prove the 'only-if part, we take the test function 

.T .T-t 

ri{t,x)— / hn{(T ~ t)ip{(T, x) da = / hn{a)Lp{a + t.x) da, t G (0, T), a; G il, (21) 



with arbitrary n G N and nonnegative (p S H2 (^t) satisfying •^\t=T = 0; 77 is nonnegative since 
Lp and hn are so (see Section [2]) . Then we have 

r]t{t,x)= / hn{cF — t)ipa{cF,x) da, a.a. (i,x) £ JIt- 
By Fubini's theorem, we have 



for all V'l, ■'/'2 e i2([0, T]). So it follows from ^ and kn = hn*k (c.p. ((H])) that 

/7 

Jo Jn 



(pt[kn * (u — uq)] + (^n * a{-,x, u,Du)\Dip) — [hn * b{-,x,u, Du)]ip] dxdt < 0, 



for all 71 G N. Observe that fc„ * (u — uq) e o-^^KPi^]; -^2(^))- Therefore, integrating by parts 
and using ip\t=T — yields 



T 

Jn 



((pdt[kn * {u — Uq)] + {hn * a(-,a;,M, I?m)|D(/j) — [hn * b{-,x,u, Du)]ip] dxdt < 0, (22) 



for all n G N and (p G i?^2' (^t) with p[t=T = 0. By means of a simple approximation argument, 
we then see that (|22|) is vahd for any ip of the form <p{t,x) = X(ti.t2)[i)^{^)^ where X(ti,t2) 
denotes the characteristic function of the time-interval (^1,^2), < ti < ^2 < T, and if) G H^iyi) 
is nonnegative. Relation (pn|) follows now from the Lebesgue differentiation theorem. D 

Our proof of the sup-bounds for subsolutions stated above relies on the subsequent truncated 
energy estimates. 

Proposition 3.1 Let p > 1, T > and fi C M^ he a hounded domain. Let the assumptions 
(K1),(K2),(Q1)-(Q5) he satisfied. Suppose that uq G L2(^) Js essentially hounded above in fl. 
Then for any weak suhsolution u G Vq,p of {I]j with esssupp^ u < 00 and any k satisfying the 
condition 

K>k := max{0,esssupuo, esssupu}, (23) 

n Ft 

there holds 

l(" - '«)+lL,([O.T];L.(n))+ \Diu - K)+lL^([o,T]xa) 

<C( f f u'dxdt+{j [A^{t)[dt)~ 

^Jo JA^it) Jo 

where 

A^{t) = {xen:u{t,x)>K}, iG(0,T), 

and the constant C = C(A^,p, q, Co, cq, C2, C2, 7, s, |^|l,([o,t]), \va + (p2\L,{nT)^T, [Q,\). 



Proof. Let u G Vq^p be a weak subsolution of JT]) in fix- Then (PO)) holds with the '<' sign 
for any nonnegative function ip e i/Hil). For t e (0,T) we choose in (pn|) the test function 
ip = u^ := (uk)+, where we set u^ '■= u — k, and k € M satisfies (P5]) . The resulting inequality 
can be written as 

(u+9t(/c„ * w„) + (/i„ * a{-,x,u,Du)\Du'^)\ dx 

i[hn*h{-,x,u,Du)]u^ + u^{uo~ Hi)kn\ dx, a.a. tG(0,T). (24) 

By positivity of fc„ and (P5)) . 

u^(uo — K)fcn rfa; < 0, a.a. i G (0, T), 



< 

In 



Thanks to (fT4ll we further have 



wj 



:9t(A:„*u,)> -at(A:„*(u+)2), a.a. (i,a:) e f^r. (25) 

Using these relations we infer from (IMl) that for a.a. t £ (0,T) 

[2 '5t[^" * (■"k)^] + (^n *a(-7 2;,u, I?w)|i:'M+)J dx < / ([hn*b{-,x,u,Du)]u'^dx. (26) 



/o ^^ ' Jn 

We next convolve (|26p with the nonnegative kernel I from assumption (Kl), and observe that in 
view of 

kn*{utfeoHl{[0,T];L,in)) 

and kr, — k * hr, we have 



I * d, 



(fcn * «)') = a,(z * fc„ * {utr) = /.„ * («+)^ 



Sending then n — > 00, and selecting an appropriate subsequence, if necessary, we thus obtain 
1 



{u'^fdx + l* / {a{-,x,u,Du)\Du+)dx <l* / b{-,x,u,Du)]u+ dx (27) 

2 Jo 7n Jn 

for a.a. i e (0,T). 

By the structure condition (Ql) we have 

(^a{t,x,u,Du)\Du^) dx ~ / (^a{t,x,u,Du)\Du) dx 

n JA^ii) 

> f (Co\Du\P - co\up - ^0) dx. (28) 

Employing (Q3) and Young's inequality we may further estimate 

/ |5(i,x,M, I?m)mJ| dx < / \C2\Du^^^u^ ^ C2\u}\'~^ul. ^ Lp2ul^\dx 

<[ &Du\P + Cs\up + ^2ut)dx, (29) 



10 



where the constant C3 > depends only on Co, C2, C2, and 7. From ([?7)) . ([25]) . and ([^ we infer 
that for a.a. t e (0, T) 

/ (u+y dx + Col* [ \Du+\P dx<2l* F, (30) 

Jn Jn 

where 

F{t) = I {{Ci + co)|ur + (/^o + (^2^+) rfa;. 



Dropping the second term in (I30p . which is nonnegative, and applying Young's inequality for 
convolutions yields 



\u 



KlL2,([0,T];L2(a)) - I("k) U,([0,r];Li(n)) 



< 2|^|l,([o,t])I-F'Ili([o,t])- (31) 

On the other hand, we may also drop the first term in pop . convolve the resulting inequality 
with fc, and use that k * I = 1, thereby obtaining 

Co|C^^+li^([o,T];L,(o))<2|^lL,([o,r]). (32) 

By Holder's inequality and assumption (Q5) we have 



/ / ^odxdt<\^o\LAnT){ [ \A^{t)\dt)'' . (33) 

Jo J A At) Jo 



Next, set 



'/3 , (l-/3)\-^_ 1-i + l 



Then the term involving ip2 can be estimated as follows, where we use (Q5), Holder's and Young's 
inequality, as well as Proposition 12. II 

ip2u^ dxdt < \ip2\L,{nT)Wt\L,.{!^T){ / \A^{t)\dt)^^ 

JA^{t) Jo 

T „ / 

+ \P u,+|i- 



•^ 



<e''|Z?u+i''^ u.+ i-^a-W 







<^''(l"«lL,([0,T];L.a.)) + l^^^;rrz.,(O.))+^""'(^''V2|l(^,)( / |A(t)|dt)^^, (35) 

J 



for every e > 0. 

Combining (PT|) -(I55 |) and choosing e such that 



26" 



(l^k([0,T]) + ^) 



1 \ 1 
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yields 

I^K lL2,([0,T];L2(n)) + \Du^ lLp([0,T];Lp(a)) 



^ / I \-^ — uy , , \u[' dxdt 

'-^0^ ^ Jo J A^(t) 



+ Wo\L.in.){ r|A.(i)|di)^+£-"'(7"V2ll(o,)( r \A,{t)\dty-^^^]. (36) 



^1 cT ^'{r-s') ■ 

\A^{t)\dt) 



Note that 

s' rs' r rj r — rj W'il^ ) 

The last condition is exactly the one in (Q5). Hence we may estimate 

(/ \A4t)\dty-^^^ <{T\n\y-^^^~^{[ \A4t)\dty, 

J{) Jo 

and thus (|36)) implies the desired energy estimate. D 

The corresponding result for supersolutions reads as follows. 

Proposition 3.2 Let p > 1, T > and il C M.^ be a bounded domain. Let the assumptions 
(K1),(K2),(Q1)-(Q5) be satisfied. Suppose that uq G i2(^) is essentially bounded below in fl. 
Then for any weak supersolution u G Vg.p of ([TJ) with essinfr^ u > — oo and any k satisfying the 
condition 

K> k := — min{0, essinf ug, essinf u}, (38) 

we have 

I (U + «)- li2,([0,T];L2(O))+ \D{u + k)_ |^^(jo.T] xn) 

<C( f f {-uydxdt+{[ \A^{t)\dt)'^) (39) 

^"'0 JA^it) Jo ' 

where 

Kit) ^{xeVL: -u(t, x) > k}, te (0, T), 

and the constant C is like in Proposition \3.1\ 

Proof. The proof is analogous to the previous one. (|20p now holds with the '>' sign and we take 
'4' = ~{u + k)- > 0. Replacing u by —u and Ai^{t) by A„(i) the same line of arguments as above 
yields the asserted estimate. D 



4 Iterative inequalities 

Let u e Vq^p be a weak subsolution of ^ in r^^. Set 



K„ = k(2- —], n = 0,1,2,..., 
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where k > niax{K, 1} will be chosen later. We further put 

Yn^ / [U- KnV+dxdt, 72 = 0,1,2, 

"'0 ^A,„„> 



'0 •^^^..(t) 

By Proposition 13. II we have for all n = 0, 1, 2 

|2 



{u - Kn+l)+lL2,([0,T];L2(fi))+ l-^(" ^ '^"+l)+lLp([o.T]xn) 



<C([ f u''dxdt+{f |^,„^^(,)|dt)='). (40) 

^Jo J A^^^^it) Jo ' 



To estimate the right-hand side, note first that 



<7 ^/ / (u-.c„rda;dt= ^—Yn. (41) 



Further, 

Yn> ( I {u~Kn)ldxdt> [ [ u'^ (l - ^^^'^ dx dt 

1 f 

> ; — —rr / U^ dx dt. 

- 2-^(^+2) L 
Hence (PO]) implies that 

(u - K„+i)+|^^^Q(, j,].^^(f2))+ |i:'(u - '««+i)+lip([o,T]xa) 



< C 2T("+2)r„ + C (^^) ' ^n ' . (42) 



On the other hand, we have by Holder's inequality and Proposition l2.1l 

Yn+l = / / {u ~ K,i+i)']_ dx dt 

<(/ {u~nn+iY+dxdty n iA,„^,(t)idi) 

~ / \/J7/ \ (1-/3)7 

< C'^(^|D(u - K„+i)+|i^([o,T]xO)j \\{u - Kn+l)+|L2,([0,T];L2(n))j 

where C = C{N,p, q) and r and f3 are given by P^ . Recall that 

!_/? (1-/3) _ |(1-|) + ^ <1 
^ P 2 l-i + ^ 
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Using (gH) and dH]) it then follows that 






kT 



>;^ 



< (cc'4^+2) \-^(i-^) (2^(^+1)) " (y„^+'"^ + ri+'"^). 

Note that (see also (I57|) ) 

7 + i > i <=? 7 > — <?4> S > — z— T-r . 

77s' r rjs' r r — rj Jt\^^ J 

The last condition is satisfied thanks to (Q5). Hence 

a ~7(— --) >0 
\r]s' r/ 

as well as 

S :=7f- - -) >a>0. 
V77 r / 

It follows from Lemma [^^ that y„ — > as n — >■ cx), provided that 

which in turn is certainly satisfied if 

/ uldxdtKC-^Ki^^-^'), (43) 

Jn 

where 

The conditions (j43p and k > maxJK, 1} are fulfilled when we set 

K := K + max|l, (C* / / u| dxdt) ^<"-"'' |. 
Since k„ — ^ 2k as n — > cx) we thus obtain 

isupu < 2k = 2(K + niax|l, (C* / / uldxdt) -''■'-'''> |j. 
Finally a short computation shows that 

^ p\ 'q) ~^ TTq 

The first part of the theorem is proved. 

The second part is proved analogously replacing u by —u and ^^(t) by Af^{t) in the previous 
arguments and employing Proposition 13.21 D 
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5 Proof of Theorem 11.21 

To prove Theorem 11.21 note first that the kernel k ~ gi-a satisfies (Kl) and (K2) with I = g^ 
and any g € (1, j^)- In particular, any q> 1 satisfying condition ((TO)) is admissible. 

Let now g > 1 be fixed such that (flO)) holds and suppose that u G Vq,p is a weak solution 
of ^ in J7t- Fix 7 S (I,??), where 77 is given by ([M)) . Note that rj < r. Hence we may apply 
Theorem 11.11 which yields boundedness of u in Qt and the estimate 



|M|Lo„(nT) ^ 2( max{|uo|L^(n),esssup|u|} + max|l,C( / / \u\'^dxdt)'' ^| 
^ Tt '^ Jo Jn ' 



(44) 



where the constant C depends only on the data. 

It remains to derive an a priori bound for the integral term involving |m| in terms of the data 
and the quantity k := max{|Mo|L^(a),esssupp^ |u|}. To this purpose we write \u\ = u+ + (^'w)+ 
and estimate |u+|l (Ut) and |(— u)4.|l (ji^) separately by means of the energy estimates from 
Propositions 13.11 and 13. 2[ respectively. 

We have 



ldxdt<2''(f f {u-K)ldxdt + K''T\n\Y 



u 
Jn '^ Jo Jn 



and by Propositions 12. II and 13. II 

'' f{u- k)1 dxdt < \iu - k)+\1^^ {TlQl)"^ 
Jn 

<e^C'^c( I f uldxdt+{T\n\)7r) +e-^{T\n\)^^^, 
for all e > 0. Choosing e sufficiently small we get a bound 

f uldxdt<C^CiN,p,q,s,\f\L^inr).T,M,f^)- 
Jn 

The bound for |(— u)_|_|x, (sir) i^ obtained analogously. Combining these estimates and (l44l) proves 
the assertion of Theorem 11.21 D 

The condition on / is sharp, at least in the cases p — 2, a ^ {0, 1), and p > 2, a G {p' /N, 1), 
as we will show in the following. 

Suppose u is a solution of equation (J9l) with smooth data mq, e.g. uq — 0. Then the optimal 
regularity for u is determined from the conditions 

d^{u - uq) G Lsiilr) and div {\Du\p-^Du) e Lsinr). (45) 

In the linear case p = 2 these lead to the maximal regularity class 

ueZ:=Hf{[0,T];L,{n))r\L,{[0,T];H^{n)), 
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where H"{J;X) denotes the vector- valued Bessel potential space of X-valued functions on the 
interval J. In fact, the first space is a consequence of the first condition, by well-known properties 
of the fractional derivation operator, see e.g. |21) . 

The question is now for which s > 1 we have Z ^^ Loo(^t)- This can be determined by 
means of cross interpolation and Sobolev embeddings. By the mixed derivative theorem (see 
[13) we have for all x? e [0, 1] 

Thus Z '^ Loo(f^T) if there exists d G (0, 1) such that 

1 N 

(1 ~'d)a> - and 2i? > — . 
s s 

This is equivalent to 

N 1 

-<^<1--. (46) 

2s sa 

There exists d E [0, 1] satisfying (|46l) if and only if s > -f- + ^, which is exactly the condition for 
s in Theorem 11.21 in the case p ~ 2. 

We now discuss a nonlinear case, namely let us assume that p > 2 and a G {p' /N, 1). Here 
we are led to the maximal regularity class 

uEZ:=Hfi[0,T];Ls{n))nLsip-i){[0,T];HUn)) with ^ = 1 + ^. 

r s N 

To see the second space suppose that u G Lq{[0, T];H?{Qj) satisfies the second condition in (HSI) . 
Assuming f < N we have DiDjU e -^^^([0, T]; Lf{^)) and 

D,u E Lq{[{),T]-Hl{n)) ^ Lq{[Q,T];L^n)), 
for all i, J = 1, . . . , n, and so the structure of the p-Laplacian leads to the condition 

ij,([0,T];L^,;«^(17)) X L,([0,r];L,(^)) C i.([0,T];i,(r!)). 

This means we need 

p-2 1 1 , {N-f){p-2) 1 1 

i-^ h - = - and -^ :^ + - = -, 

q q s rN r s 

which in turn implies q = s{p — 1) and f = r. 

As before we have to determine those s > 1 for which we have Z ^^ Loo{^t)- Note first that 
the conditions a > - and 2 > -^ are necessary. Assuming this it follows that s > max{ — , — }, 
in particular we have s > f due to p > 2. By the mixed derivative theorem we then have for all 

^ e [0, 1] 

Z -> i?,(p:i)^<''-" ([0,T];i,(r!)) n L,(p_i)([0,T];iJ. ^ + = (O)) 
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with sharp embeddings. Hence Z ^^ Loo{^t) if there exists z9 G (0, 1) such that 

_,,/! 1\ 1 / N N\ N 

{l-^)[a~- + - -T >^ -T and z? 2--r + - >-, 

V s s[p—l)/ s(p — 1) \ r sJ s 

which is equivalent to 



r s s s{'p—l) 

Thus it boils down to the condition wi < aj2- A short computation shows that this condition is 
in fact equivalent to s > — + ^ . 

Recall that we assumed that f < N, that is f = f < A^. The condition f < A^ in turn is 
equivalent to s < A^. Thus s has to satisfy the condition 

N 1 

h - < s < iV, 

p a 

which is possible, by the assumption a e [p'/N,l). Hence in this case the condition on / is 
optimal as well. 

6 Homogenous structures 

In this section we consider the special case of homogenous structures. By this we mean equations 
of the type 

dt{k*{u-uo)\ -d:iYa{t,x,u,Du) = {), te(0,T),xeO, (47) 

where 

(HS) (a(t,a;,C,?7)|r/) > C^H^ , Ht,x,^,7^)\ < CM"-', 

for a. a. {t,x) G Ht, and all ^ G M, 77 G M^. Here Co and Ci are positive constants. 

In this situation the weak maximum principle takes the same form as in the classical parabolic 
case. Moreover the assumption (K2) can be dropped, g = 1 is here admissible. 

Theorem 6.1 Let p > 1, T > 0, and fl C M^ be a bounded domain. Assume (Kl) and (HS), 
and let uq G L2(^)- Then for any weak subsolution (supersolution) u G Vi,p of |^7p , we have for 
a. a. (t, x) G ^T 

u{t, x) < max-^ O,esssupwojesssupu > { u{t,x) > min < 0,essinf uojessinf u > 

provided this maximum (minimum) is finite. 

Proof. Note first that Lemma 13.11 also holds under the assumptions of Theorem 16.11 It suffices 
to consider the subsolution case. We take 

K = k = max<^ 0, ess sup Mo, ess sup u >, (48) 
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assuming that this quantity is finite, and proceed as in the proof of Proposition 13. II This yields 

l(w - ^)+lL(o.) + l^(" - ^)+li(o.) < 0, 
which immediately implies the assertion. D 



Theorem lG.ll shows in particular that the maximum principle holds in the usual form for weak 
solutions of the time fractional p-Laplace equation ^ with / = 0. 

We remark that the case q — I can occur. In |20[ Section 3] an example is given for a kernel 
k satisfying (Kl) with / ^ Lq{[0,T]) for all q > 1 and T > 0. 

7 Natural growth conditions 

Finally we consider the case of 'natural' or Hadamard growth conditions with respect to \Du\. 
For the sake of simplicity we suppose that 

(Q) ia{t,x,^,TiM>Co\r,\P, \ait,x,^,i^)\ < Ci\r,\P-\ \b{t,x,^,r,)\ < C2\v\p, 

for a. a. {t,x) G fir, and all ^ G M, 77 G M^, where Ci, i = 0,1,2 are positive constants. In the 
classical parabolic case it is known that weak solutions of the corresponding problem under the 
conditions (Q) are in general not bounded. However there are results (also in a more general 
situation) which provide Loo bounds in terms of the data assuming in addition that the weak 
solution is bounded, see e.g. [7]. It turns out that corresponding results can be obtained for ([T]). 
Here we only prove such a result in the case where (Q) holds. It generalizes Theorem 4.3 in [20] . 
where p = 2 is required. As in the previous section we may drop assumption (K2). 

Theorem 7.1 Let p > 1, T > 0, and fl C M^ be a bounded domain. Let (Kl) and (Q) be 
satisfied, and suppose that uq G Loo{^)- Then for any bounded weak solution u G Vi^p of (Qp, 



\u\L^{nT} < max||uo|L„(o),esssup|u|l 



Proof We proceed as in the proof of [3 Theorem 17.1], see also [20, Theorem 4.3]. Assume that 
K :— esssupQ^ u > k, where k is as in ([^S]). Let e > be such that k := K — e > k. We then 
choose the test functions w^ = (u — k)+ and estimate similarly as above, using the conditions 
(Q). This yields 



k;tlL(n.) + l^^«li,(o.) < CiCo,C2) \Dutfut 

<eC{Co,C2) \Du+\" 



Li(nr) 



Li(Ot) 

Choosing e sufficiently small, it follows that |uj|| (o ■) < 0, that is u < k < K a.e. in ftx, a 
contradiction. Hence u < k a.e. in J7y. The lower bound is obtained analogously. D 
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